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Abstract 

We investigate the thermal quenching of the multimodal photoluminescence from InAs/InP (001) 
self-assembled quantum dots. The temperature evolution of the photoluminescence spectra of two 
samples is followed from 10 K to 300 K. We develop a coupled rate-equation model that includes 
the effect of carrier thermal escape from a quantum dot to the wetting layer and to the InP 
matrix, followed by transport, recapture or non-radiative recombination. Our model reproduces 
the temperature dependence of the emission of each family of quantum dots with a single set of 
parameters. We find that the main escape mechanism of the carriers confined in the quantum dots 
is through thermal emission to the wetting layer. The activation energy for this process is found 
to be close to one-half the energy difference between that of a given family of quantum dots and 
that of the wetting layer as measured by photoluminescence excitation experiments. This indicates 
that electron and holes exit the InAs quantum dots as correlated pairs. 

PACS numbers: 
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I. INTRODUCTION 



Quantum dots (QDs) are designable mesoscopic atomic assemblies whose effective 
electronic density of states is 5-functionlike.- One of the most studied systems is self- 
assembled QDs grown in the Stransky-Krastanov mode. It has been demonstrated that 
self-assembled QDs can find applications in fields ranging from nano-optoelectronics^ to 
quantum computing. 3 Understanding the processes that result in the thermal quenching of 
the photoluminescence (PL) of QDs is thus important not only on fundamental grounds but 
also for the realization of efficient devices operating at room temperature. 

In the case of InAs/GaAs QDs, it is now well established that the two mechanisms that 
control the populations of electron-hole pairs in the ground state of QDs are their radia- 
tive recombination and their escape to higher lying energy states. Thermal quenching then 
results from non-radiative recombination processes that occur in one or several of those 
higher energy states.-^^&^ii^ 1 ^ Even though both measurements and theoretical mod- 
eling appears straightforward, there remains to this day significant differences and apparent 
contradictions in the interpretation of the results published by different groups in the last 
decade. 

Two important questions remain unanswered. First, the identification of the higher 
energy states that contain non-radiative centers. Potential candidates are (i) the so-called 
wetting layer (WL), which is a few monolayer-thick InAs pseudomorphic quantum well (QW) 
always present in the Stranski-Krastanov growth mode; (ii) the matrix in which the QDs 
are inserted, either GaAs or a confining QW; (iii) defects or impurities in the matrix; or 
(iv) QD excited states. Second, it is still not clear whether the confined electron and holes 
escape a QD as a unit (exciton), as a correlated e-h pair or as uncorrelated electrons and 
holes. The e-h correlation can be evidenced by the values of the activation energies 

Ef = isAEi , (1) 

where AEi is the difference between the energy of the higher energy state i and that of the 
QDs. If v = 1, e-h pairs escape as excitons; if v = 1/2, the escape mechanism involves 
correlated e-h pairs, while if v < 1/2, it involves uncorrelated electrons and holes.- 

Experimentally, QD-size fluctuations result in a distribution of their quantized ener- 
gies and hence to a distribution of activation energies. As temperature is increased, re- 
distribution of the carrier population occurs towards the low energy tail of the QD energy 
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distribution.-^^ When the full-width at half-maximum (FWHM) of the emission of the 
ensemble of QDs is much smaller than any expected activation energy, the QD energy dis- 
tribution is often represented by a 5-function. The solution of the thermal rate equations in 
steady-state for the ensemble photoluminescence (PL) intensity / then gives 



where m is the number of high energy states involved and is the activation energy for 
the transfer from QD to state i.— Using this approach, Le Ru et al. identified the GaAs 
matrix as the dominant recombination state and obtained v = 1.- Using samples where the 
QDs were imbedded in an InGaAs QW, Torchysnka et al. found several activation energies 
that corressponded well to transitions from QDs to QW and from QW to the GaAs ma- 
trix, assuming v = 1.— ^ Seravalli et al. examined samples where the QDs were inserted 
in In^Gai-^As and In^AlyAs confining layers. They found two dominant activation ener- 
gies that pointed to exciton transitions from QDs to an unidentified low-energy defect and 
transitions from QDs to the WL.— 

More information is available if the internal population redistribution is also measured 
and modeled. In that case, the rate equations must take into account the QD energy 
distribution. Yang et al. analysed the PL decay times within the QD emission band with a 
model that only included the WL as the high energy state. 4 They found a good agreement 
with experimental data with v ~ 1/2. Sanguinetti et al. developed a model that included 
exciton transitions from QDs to WL and WL to the GaAs matrix.- It reproduced well the 
QD PL integrated intensity, FWHM and peak energy as a function of temperature. Dawson 
et al. included in their model independant distributions of electrons and holes but they only 
considered transitions from QDs to the GaAs matrix.- 1 ^ They found that the PL integrated 
intensity, FWHM and peak energy was best fitted by uncorrelated carrier escape. 

The reason for these discrepancies is unclear. One might point out to different escape 
mechanisms in different samples, to insufficient information provided by a single QD emis- 
sion band, or to differences in the theoretical models and/or assumptions underlying their 
treatment. On the other hand, little work of a similar nature has been carried out on equiva- 
lent QDs grown on materials systems with an energy spectrum quantitatively different from 
that of InAs/GaAs. In order to solve some of the problems mentioned above, we present a 
study of the thermal quenching from two samples that contain InAs QDs embedded in an 
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InP matrix. Their PL spectra show several well resolved emission bands, extending from 
0.75 to 1.1 eV, that can be associated with families of QDs that differ in height by one mono- 
layer (ML).— ii£ii&iLi&i£i22j2I The evolution with temperature of this multimodal PL imposes 
stringent constraints on a model based on coupled rate equations as it should reproduce the 
thermal behavior of many peaks with the same set of parameters. 



II. EXPERIMENTAL PROCEDURES 



Sample A was grown by low-pressure metalorganic vapor phase epitaxy on (001) InP Fe- 
doped semi- insulating substrates at a reactor pressure of 160 torr. After the growth at 600 
°C of a 100-nm InP buffer layer, the temperature was lowered to 500 °C during 90 s while 
still growing InP. 2.4 InAs MLs, sandwiched between two 20-ML thick InP layers, were then 
deposited at 500 °C. The temperature was raised back to 600 °C and the growth terminated 
by the deposition of a 30-nm InP cap layer. Trimethylindium (TMI), tertiarybutylarsine 
(TBA), and tertiarybutylphosphine (TBP) were used as precursors, and Pd-purified H2 as 
the carrier gas. The flow rates of TMI and TBA were kept at 0.05 and 0.95 /imol s _1 and 
the TBP flow rate was 3.5 /zmol s _1 . The growth rate of the InAs layers was close to 0.4 
ML s -1 (0.38 fim h _1 ). The growth interruption sequence used to switch from InP to InAs 
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In particular, a 4s growth interruption was applied after 



and back is described in Ref[: 
InAs deposition. 

Sample B was grown on (001) InP substrates by chemical beam epitaxy from TMI, arsine 
(AsH 3 ), and phosphine (PH 3 ). AsH 3 and PH 3 were cracked at 850 C in a fast switching 
high-temperature injector to produce predominantly As2 and P2. After desorption of the 
surface oxide, the growth was initiated with an InP buffer layer. Then, TMIn and As2 were 
injected simultaneously into the chamber to grow about 2.2 MLs of InAs layer. This was 
followed by a 30 s growth interruption time under an overpressure of As2 to allow the QDs 
to form. The QDs were then capped with InP. Further information on the growth procedure 
can be found in Ref 21 



The PL measurements were carried out with the samples mounted strain free in a helium- 
flow cryostat. The excitation source was the 632.8 nm line of a He-Ne laser. The signal was 
analyzed with a spectrally calibrated, nitrogen-purged Fourier transform spectrometer us- 
ing a liquid-nitrogen cooled InSb detector. The PLE was excited with the monochromatized 
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light of a 150 W tungsten lamp. The PL was then analyzed with a 0.5-m double grating 
spectrometer and detected with a liquid-nitrogen cooled InGaAs photodetector array sensi- 
tive up to 2 /im. The PLE spectra were not corrected by the wavelength dependence of the 
excitation intensity. 

III. RESULTS 

Figure Q] shows the low-temperature PL and PLE spectra of the samples. The emission 
of sample A comprises five peaks while that of sample B comprises nine peaks. The energy 
position of peak Bl* encompasses that of peaks Al and A2 while that of peaks B3 to B5 
are close to that of peaks A3 to A5. 

The PLE spectra of peak Al shows an edge at 1.19 ± .01 eV, labeled WLhh, that is close 
to the exciton energy in thin InAs QWs.— ^ It can thus be attributed to the WL. It follows 
that the high energy tail of the emission from sample A at low-temperature corresponds to 
residual emission from the WL. 

The PLE spectra of peaks B3 and B6 are also shown in fig. [TJ The low-energy edge in 
both spectra can be attributed to the first excited state QD^. Another edge appears in both 
spectra at the same energy of 1.30 ± .01 eV. Actually, this feature is common the the PLE 
spectra of all peaks from sample B. It can thus be associated with the WL^ transition in 
sample B. We attribute its higher energy with respect to that of the WL in sample A to the 
longer interruption that took place during the growth of sample B. This probably allowed 
the formation of thicker QDs and hence, a thinner WL with a blue-shifted resonance energy. 
The difference in energy of the fundamental WL optical transition is useful for our purpose 
as it adds another constraint to the thermal model described in Section IIVI 

The evolution of the PL intensity of both samples as a function of temperature is depicted 
in Fig. [2j The emission from peak Al is rapidly quenched for T > 100 K while that of 
peak A2 remains nearly constant for T < 170 K. The intensity of peaks A3 to A5 actually 
increases for T ~ 200 K before decreasing at higher temperatures. The emission from sample 
B is more robust as only peaks Bl* and B3 show a significant intensity decrease at 300 K. 
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FIG. 1: Low-temperature PL and PLE spectra of samples A and B. QDih, W^, W/^, and InP 
refer to excitonic heavy- and light-hole resonances in QDs, WL and the InP matrix respectively. 



IV. RATE-EQUATION MODEL 

Our thermal model is similar to those developed in Refs |4j and [5|. A series of coupled 
steady-state rate equations that control the populations rii in each state i is obtained from 
the detailed balance principle:- 



- Thi 



rijUji + Pi = , (3a) 

where Ri is the recombination rate of state i, Ni the number of states per unit area, and Uij 
the transfer cross-section from state i to state j. Pi represents the carrier generation, is 
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FIG. 2: Evolution of the PL spectra of samples A and B as a function of temperature, 
defined in Eq. ([Q) and 

^ = ^exp|^ ( ^ fc ~ Ej) | if Bk >£,-. (3b) 

We improved on previous models by incorporating several unique features. 
(i) As both WL and InP matrix features are observed in the PLE spectra, transfer from 
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QDs to WL, from QDs to the InP matrix, and from WL to the InP matrix are allowed. 
(ii) The allowed states of the WL and the InP matrix are distributed over a wide range of 
energy above their fundamental edge, for the WL and E^ 1 for InP. To make this fact 
numerically tractable, the WL is separated into segments of width AEw extending from 



to Eg' 1 . The effective number of states per unit area of each segment is N\y 
where Dw 



D]yAE\y , 



,InAs 



m x '"^ s /nh 2 is the two-dimensional density of states in the WL. The InP matrix 



is similarly segmented. We assume a ■sjE — E g dependence for its three-dimensional density 
of states. The number of states per unit area for a segment that extends from E™ m to 

E min + ig thug 



(4) 



N Mi = -D M £ M [(Er n + AE M - Ef f 2 

where D M = 2 1 ^ 2 (m I x nP ) 3 ^ 2 /(ir 2 h 3 ) (Ref. [2^) and £ M is the active thickness of the matrix. 
(Hi) E^ (T) and E^(T) are assumed to follow the Varshni temperature dependence with 
the parameters of bulk InAs and InP, respectively.— 

In order to restrain the number of adjustable parameters, the following assumptions were 
made. 

R D if % e D, 

Rw for W lowest energy segment, 
Rm for M lowest energy segment, 
for all other states, 



Ra 
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P; 



u D if i e D, 
1 if i e W, 

P if i = M highest energy segment, 
if not, 
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(5d) 



where D, W, and M refer respectively to the ensemble of QD, WL, and InP matrix states. 
We further assumed that the parameters in Eqs [5] are independant of temperature. 

In Eq. ISU Re> corresponds to QD radiative rate. As no emission from the WL nor 
the InP matrix is observed at high temperature, R\y and Rm correspond to non-radiative 
rates. R\y was assigned only to the lowest energy segment because, in a QW, free excitons 
form 2-dimensional polaritons that do not couple to photon-like polaritons propagating 
perpendicular to the QW plane unless their energy is within a small bandwidth near the 
bottom of the band.— Excitons must thus relax to the bottom of their energy band before 
they can recombine radiatively or non-radiatively. 27 A similar argument can be made for 
Rm-— Figure [3] schematizes the rate-equation model used to analyze our data. 
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FIG. 3: Schematic of the coupled rate-equation model. 
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TABLE I: Parameters obtained by fitting the model to the experimental data. The numbers in 
parenthesis are the uncertainties calculated with Eq. [6J 





Rw/Rd 


Rm/Rd 


hi 


UwdDq/Rd 


UmdDq/Rd 


U MW vd £f(0) 








(nm) 






(Hz cm 2 ) (eV) 


Sample A 


22 (1) 


9 (18) 


100 a 


12.3 (0.5) 


5000 (300) 


0.71 (0.03) 0.64 (0.03) 1.19 (0.04) 


Sample B 


36 (2) 


20 (10) 


100 a 


51 (2) 


20 (140) 


10 (10) 0.51 (1) 1.30 tt 



"Fixed parameter. 



V. DISCUSSION 

To compare the simulations with the experiments, the data were treated as follows. 
The peak energy and integrated intensity of each of each peak at a given temperature was 
obtained by fitting the PL spectrum with a series of gaussians peaks. This procedure was 
found to reproduce well the PL spectra except for peak Bl*. Its peak energy and intensity 
was obtained by substracting the intensity of all the other peaks from the total intensity of 
the PL emission. The energy position of each family of QDs served as input to the model. 

In the model, the total number of QD states per unit area Do and the QD recombination 
rate Rp are scaling factors. The relative number of states for each QD family was assumed 
to be given by the relative intensity of the PL at low temperature. AEw and AEm were 
set at 10 meV, a value close to the spectral extent of the absorption edge. In Eq. l5dl the 
parameters Uww and Umm were set to a high value to ensure that the excitons in the WL 
and the InP matrix are in thermal equilibrium. Finally, we fixed Im at 100 nm, a value 
close to the penetration length of the excitation source.— The model is thus left with seven 
adjustable materials parameters: Rw/Rd, Rm/Rd, UwdDq/ Rd, UaidDq/Rd, Umw, u d, 
and Ef (0). 

The result of our simulations is presented in Fig. HI The temperature dependence of all 
peaks from sample A are very well reproduced by our model with the optimized parameters 
listed in Table [H In particular, our model reproduces the intensity increase of peaks A3 to 
A5 when T > 180 K. The uncertainties Aa^ of the optimized parameters were estimated 
using^ 




Stricktly speaking, Eq. [6] is only valid if the cross-partial derivatives dx 2 /daidj are small 
with respect to the diagonal terms, which is not the case here. However, it gives a good 
estimate of the sensitivity of the fit to a given parameter. 

It can be seen from Table [I] that the parameters with high uncertainties are relative to 
the InP matrix, an indication that the main QD escape channel is through the wetting layer. 
To further test this hypothesis, we deactivated the contribution to the thermal quenching of 
InP matrix by fixing Rm = 0, with no significant change to the fit. On the other hand, no fit 
could be achieved when the the WL was similarly deactivated. Further, the fitted value for 
E^(0) corresponds within uncertainties to the measured value of the low energy edge WLhh 
shown in Fig. [TJ We can thus conclude that in sample A, the main quenching mechanism 
is through carrier escape from QD to WL followed by a non-radiative recombination of the 
carriers in the WL. 

There is globally much less thermal quenching in sample B and thus less dynamics to 
constrain the model. To extract relevant information, we fixed £^(0) to the value of WLhh 
obtained the PLE spectra of sample B, as shown in Fig. [U Here also, peaks Bl*, B3, and 
B4 are well reproduced by our simulations. There was nearly no change of intensity for 
peaks B5 to B9, while our model predicts a slight increase. The discrepancy can easily be 
explained by our neglect of the temperature dependence of the parameters of the model. We 
also note that the two most significant materials parameters, Ry/ and Uwd have realistic 
and comparable values for both samples. 

As for the parameter ud, the simulations indicate that it is close to one-half for our 
samples. Further, our model could not reproduce the data if vn was fixed to 1. This 
indicates that electron and holes escape from the QDs mostly as correlated e-h pairs. We 
thus corroborate the findings of Yang et ai— 

It is instructive to simulate with our model the temperature dependence of the integrated 
intensity of a monomodal QD emission. We have generated a gaussian distribution of fifteen 
QD subfamilies centered at an energy E D and shifted with respect to the WL by AE = 
Eg —E D . The FWHM of the distribution was fixed at 0.25A.E, a value typical of monomodal 
InAs/GaAS QD emission.- We used in the simulations the same materials parameters as 
those found for sample A except for ujj = 0.5. 

The result of our simulations for AE = 100, 200, and 300 meV is shown as symbols in 
Fig. The curves were analyzed with a sum of activated processes as described by Eq. 
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[2j All curves are well reproduced with the activation energies given in Table [TTJ The rapid 
drop of intensity at elevated temperatures is controlled by Ef. In all three cases, we find 
Ef > uAE. Further, Ef does not correspond to any difference between the energy levels 
present in the model. 



TABLE II: Parameters obtained by fitting the simulated data of Fig. [5] with Eq. [2j 



AE 


Ef 


E% 


(meV) 


(meV) 


(meV) 


100 


160 


50 


200 a 


190 




300 a 


230 





"Only one activation energy required. 

These simulations show that the difference between Ef and the activation energies in- 
serted in the model comes from carrier transport in the WL and recapture by QDs. This 
induces a redistribution of the carriers within the subfamilies that slows down thermal 
quenching, resulting in an ensemble effective activation energy higher than actual ones. 
Therefore, in systems where recapture competes with recombination, Eq. [2] gives empirical 
activation energies that might not correspond to any physical process at play. 

VI. CONCLUSIONS 

We have developed a system of coupled rate equations for the temperature dependence of 
the multimodal PL of InAs/InP QDs. The model includes carrier escape to the InAs wetting 
layer and to the surrounding InP matrix as well as carrier transport and retrapping. Even 
though our model comprises seven adjustable parameters, the constraints imposed by the 
simulation of the complex temperature behavior of up to five different QD families makes our 
fits robust. We find that the main quenching mechanism is induced by carrier escape to the 
wetting layer followed by non-radiative recombination. Further, our results clearly establish 
that, for both samples examined, electrons and holes are emitted as correlated pairs rather 
than excitons. Finally, we show that carrier redistribution within the QD energy levels as 
temperature in increased can yield activation energies obtained from analyzing PL integrated 
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intensities that do not correspond to any actual physical process. 

We cannot assert whether correlated-pair escape is characteristic of self-assembled QDs 
or specific to our samples. The latter case could mean that the temperature dependence of 
QD optical emission is governed by microscopic parameters such as the size and shape of 
individual QDs. A better theoretical understanding of the interactions between QDs and 
their environment is thus not only of great fundamental interest, but could also impact the 
design of QD-based devices. 
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FIG. 5: Arrhenius plot of the temperature dependence of the integrated intensity of simulated 
monomodal QD emission for three values of AE = E^ — E D . The solid lines are best fits using 
Eq. m 
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